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Abstract. We made a detailed study of the properties of the electron layer near the quark surface of strange 
stars with strong (~ 10 14 — 10 17 G) magnetic helds. The electrostatic potential and the electric held at the quark 
surface were calculated as functions of the magnetic held intensity of bare strange stars. Using an ultrastrong 
(B > 2.5 x 10 1B G) magnetic field, we found that the distribution of electrons becomes an exponential function of 
radial distance, which is quite different in a magnetic field-free case. We also calculated the variation in gap width 
between the strange core and the normal nuclear crust for strange stars, which is due to magnetic field effect. 
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1. Introduction 

Strange quark matter (SQM), made up of roughly equal 
numbers of up, down, and strange quarks, may be the ab- 
solute ground state of strong interaction. If this hypothesis 
is correct, strange quark stars (strange star, SS) may ex- 
ist. Of course, it is very important to identify SSs. A new 
method for distinguishing SSs from neutron stars has been 
proposed recently ba sed on the strong electric field at the 
qu ark surface of SSs tXu et all kOOlHlIsovLboOlah . 

lAlcock et all l)l986|) first pointed out that a very strong 
(~ 10 17 V/cm) electric field may exist near the bare quark 
surface. Due to the strong electric field, pairs could 
be created spontaneously in a few empty quantum states; 
these pairs, which then are thermalized in the electron 
layer subsequently annihilate into photons l|Usovl Il998l 
2001b). This mechanism could play an important role in 
producing the thermal emission of bare SSs, because the 
SQM is a poor radiator of thermal photons at frequencies 
less t han its plasma fre quency (~ 20MeV)dAlcock et a " 
1986). Previous studies ( Alcock et alll986t|Kettner et a 



199,'iHHuang fc LuUl997 ) have also shown that the strong 



electric field may be decisive for the existence of a nuclear 
crust on the quark surface of SSs, which may be created in 
a supernova or by accre t ion fr o m the interstellar m edium 
l|Glendenning fc Webeii Il992t Idendenning et at Il995|) . 
This is due to the fact that the strong positive Coulomb 
barrier prevents atomic nuclei bound in the nuclear crust 
from coming into direct contact with the strange core. A 
SS with a nuclear crust may have different cooling time 
scales and surface properties from a bare SS. Thus the 
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study of the nuclear crust in a SS is important for future 
observations of strange pulsars. 



As pulsars, SSs are strongly magnetized. From a 
sample of more than 400 pulsars, the surface mag- 
netic field intensity lies in the in terval ~ 2 x 
10 10 -10 13 G ^Manchester fc Tavorlll98ll). Several a uthors 
llBisnovatvi-Kogarl Il993t [Duncan fc ThompsonL Il992t 
IThompson fc Duncanl Il993h have proposed two different 
physical mechanisms that lead to an amplification of some 
initial magnetic fields in a collapsing star, fields as strong 
as B ~ 10 14 — 10 16 G, or even more. Considering the flux 
throughout the stars as a simple trapped primordial flux, 
the internal magnetic field may go up to B ~ 10 18 — 10 20 G. 
The energy of a charged particle changes significantly in 
the quantum limit if the magnetic field is equal to or 
greater than some critical value B c = rnfc 3 / (qih)G, where 
m% and qi are the mass and charge, respectively. For u 
and d quarks of current mass 5 MeV, B c ~ 4.4 x 10 15 G, 
while for s quark of current mass 150MeV, it is ~ 10 18 G. 
Especially for electrons of mass 0.5MeV, a strong mag- 
netic field that exceeds the critical value B c ~ 4.4 x 10 13 G 
forces the electron to behave as a one-dimensional gas 
rather than three-dimensional. Therefore, the quantum- 
mechanical effect of the magnetic field on the distribution 
of electrons near the quark surface of SSs cannot be ne- 
glected, just as the temperature effect leads to a consider- 
able red uction of the electro st atic potential a t the quark 
surface ijKettner et a | 119951) . IPhukonl |2000) performed 
a perturbative expansion of the thermodynamics to the 
lowest order term in B 2 / ' [x 2 to reveal the magnetic effect. 
However, in his work, the details are lost in the middle 
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magnetic field intensity, and the approach is invalid for 
the ultra-strong field case. 

In this work we present a detailed investigation of the 
properties of the electron layer in the presence of a mag- 
netic field and probe the magnetic field dependence of the 
gap width between the crust and the quark core. In the 
next section we describe the thermodynamic properties 
of bulk SQM in the presence of a strong magnetic field 
at zero temperature. In Sect. 3 we extend the Poisson 
equation near the quark surface to include the effect of 
magnetic field, and determine the electrostatic potential 
numerically. For an ultra-strong magnetic field, we sim- 
plify the equations and solve them analytically. In Sect. 4 
we present the results of SSs with a thin crust. In the last 
section, we give our conclusions. 

2. Bulk SQM in strong magnetic fields 

For a constant magnetic field B, the single particle(species 
i) energy eigenvalue is written as 



PI 



(2n + s + l) qi B, 



(1) 



where n =0,1,2,..., are the principle quantum numbers for 
allowed Laudau levels, s = ±1 refers to spin up(+) and 
down(-), and pi and qi are the components of particle mo- 
mentum along the field direction and the absolute value of 
the electric charge, respectively. Setting 2n + s + 1 = 2i/, 
where v =0,1,2,..., we can rewrite the single particle en- 
ergy eigenvalue in the following form 



£,, = 



P< 



2vq l B. 



(2) 



Now it is very easy to show that the v = state is 
singly degenerate, while the others with v ^ are dou- 
bly degenerate. To a zero temperature approximation, the 
thermodynamic potential o f the ith species is given by 
llChakrabartv fc Sahul [l996) 



if b Vi0 {v-i\fy?i - mf - 2v l q l B 

— n ^ 



(m 2 + 2u iqi B)ln 



-2viqiB 



where and are the degeneracy (equaling 6 for quarks 
and 2 for electron) and chemical potential of the ith 
species, fe^o = 1 — \8 Ui Q. The upper limit v i m is the great- 
est integer not exceeding (fif — m 2 )/2qiB. 

For a SS, the net positive charge of the quarks will 
be balanced locally by electrons up to radial distances 
T < R m , where R m is only slightly smaller than the stellar 
radius R. This charge neutrality condition of bulk matter 
is expressed as 

2n u - n d — n s - 3n e = 0. (4) 
And the baryon number density of the system is given by 

O, (5) 



1 



which is considered as a constant parameter. For bulk 
SQM, assuming the condition of /3 equilibrium, we have 



Md = Ms = M, 
M« = M - Me- 



(6) 
(7) 



To solve the above equations for a given magnetic field 
and baryon number density, we need the relationship be- 
tween rii and /ij. From the well-known thermodynamic 
relation, the expression for the nu mber density of the ith 
speci es (i=u,d and e) is given by ijChakrabartv fc SahuL 
1996) 



9iq%B 
2tt 2 



KoJtf-ml-2 mi B. (8) 



To the weak magnetic field limit (/if — mj)/2qiB — > oo, 
taking — > J dvi , we can obtain the classical expres- 

Vi=0 

sion of the number density 



*i(M? 



f)3/2 



6tt 2 

which is valid for s quark due to its large mass 



^) 3/2 



(9) 



(10) 



On the other hand, for an ultra-strong field B that satisfies 
{jif — mf)/2qiB < 1, only the state Vi = is allowed. We 
then can get rid of the sum for Laudau levels in Eq. (8) to 
read 



9iqjB 
4tt 2 



(11) 



Now we can solve Eqs.(4-7) to obtain the number densities 
and the chemical potential of all species in the interior 
including the vicinity of R m (r < R m ) of SSs with various 
values of B. 

3 Bare quark surface 



(3) 3.1 Equations and numerical results 



Beyond R m , in the region R m < r < oo, the charge neu- 
trality is a global rather than a local conditi on. A concise 
model for this layer was first developed by lAlcock et at 
(1986). In a simple Thomas-Fcrmi model, the electrostatic 
potential near the quark surface is described by Poisson's 
equation 



d 2 V 
dr 2 



4:Tve 2 (n e 
47re 2 n e , 



i g ), r<R, 
r > R, 



(12) 



where V/e is the electrostatic potential. The chemical 
equilibrium of the electrons in the layer implies that the 
value /loo = Me — V is constant. Since far outside the star, 
both V and /i e tend to 0, it follows that //qo = and 
/i e = V, and 



n q = -n h 



1 



-n d 



1 



(13) 
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is the quark charge density inside the quark surface. 
Because R and R m differ only by a few hundred fermi, 
density n q in this range can be treated as being indepen- 
dent of r. Because of the local electric charge neutrality 
inside R m , 
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n e {Rm) 



9 e eB 
2tt 2 



u=0 



>lA) 



2veB, 



(14) 



> 

CD 



where V q is clearly the chemical potential corresponding 
to the net charge density of quarks inside R rn , which can 
be obtained from Eqs.(4-7). 

Assuming the magnetic field in the electron layer is 
uniform and substituting Eq.(8) into Eq.(12), we come to 
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V(R) ; f~ ^""^ 
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d 2 V 
dr 2 



E KaVv 7 

v=0 



2veB 



B/G 



= < 



- E * 



v 



V 2 



E K^V 2 -m 2 -2veB, r>R 



2 — 2v' eB I r < R, (15^ig- 1- V q and V(R) vs. magnetic field intensity B for 

riB — 2.5no. The dotted curve shows the number of the 
maximum Laudau level v m with different B. 



v=0 



The boundary conditions for the above equation are V — > 
V q and dV/dr as r -> 0, and V -► and dV/dr -> 
oc, respectively. By integrating Eq.(15) we can 



as r - 
obtain 



dV 
dr 



-J 2 ^ E b v0 V^V 2 -ml- 2veB 



-V„ 



V 2 



2veB 



> 

a) 



-(ml + 2veB) In 



-2ueB 



-2 E h' (v 

v'=0 



v q ) 



V 2 



i 2 - 2v'eB 



r<R, 



(16) 




2M p- \ E K 

(m\ + 2veB) In 



-500 



V^V 2 



2veB 



500 1000 

Radial distance / fm 



1500 



V+y/V 2 -rni-2ueB 



r>R. 



The global charge neut rality determines the potential 
at the quark surface to be ijKettner et all Il995|) 



V(R) = V q 



P e (i? m )-P e (oo) 



Fig. 2. Electrostatic potential of electrons in the close 
vicinity inside and outside of the quark surface of Bare 
SSs. The location of the surface is indicated by the ver- 
tical line. The figures assigned to these curves refer to 
magnetic field intensity. 



(17) 



where P e (R) is the kinetic pressure of electrons, which 
could be obtained by the thermodynamic relation P e = 
— 57 e . In Fig.l we show the variation of V q and V(R) with 
the magnetic field intensity B. We can see V q is nearly 
independent of B for our range of fields (~ 10 14 — 10 17 G). 
And the relation V(R) ~ jV q applied widely in previ- 
ous studies is a good approximation for a weak (B < 
10 15 G) magnetic field, while V(R) ~ \V q for an ultra- 
strong (B > 2.5 x 10 16 G) field. In the middle range 
(B - 10 15 - 2.5 x 10 16 G), V{R) oscillates with magnetic 
field intensity. These oscillations are easily understood in 



accordance with Laudau levels. As is known, the maxi- 
mum level f m , the integer of the ratio ^"^Jb^ ' takes sep- 
arate values as a function of B. We also show the variation 
of v m with B in Fig.l. It is obvious that the oscillation 
of V(R) occurs synchronously following the skip of v m . 
For a weak field, v m has a very large value as — ► so 
that the oscillation amplitude is very small. However, the 
quantum effect becomes important due to the low value 
of v m that has the comparable magnitude of Av in the 
middle field intensity. Therefore V(R) oscillates consid- 
erably in this range depending on the separate values of 
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B/G 



Fig. 3. Electric field at the quark surface vs. magnetic 
field intensity B. 



v m . The oscillation disappears until v m vanishes for very 
strong magnetic fields. Furthermore, not only V(R), but 
also other quantities involving the summation of Laudau 
levels show oscillations (see subsequent figures). 

The behaviors of V(r) that vary with respect to B are 
exhibited in Fig. 2. We can see the stronger the magnetic 
field is, the larger the potential changes with r. To illus- 
trate the details, we plot the curve of the electric field 
(at the quark surface E(R) = {—dV/dr) r= n) versus mag- 
netic field in Fig. 3. With the increase in B, the electric 
field increases remarkably after oscillations within a range 
(~ 10 15 - 2.5 x 10 16 G). 

3.2 Ultra-strong magnetic field approximation 

At ultra-strong magnetic field (B > 2.5 x 10 16 G for 
V q ~ 17MeV), only the lowest Laudau level is allowed, 
and Eq.(16) becomes to 



dV 
dr 



2e 3 B 



V^V 2 - m\ - (2V - V q )^V 2 - ml 

r<R, 



-mi m 1 * 



2e 3 B . 



(18) 



V \/V 2 — m 2 — m 2 In 



r>R. 



Since V 3> m e , the above equation can be written more 
simply, if we ignore the mass of the electron, to read 



By integrating this equation and using the relation 
V(R) ~ \V q (obtained from Eq.(17) for ultra-strong 
field approximation), we get the expression of V, and 
E = —dV/dr as a function of radial distance r 



V(r) 



E{r) 



V q \ 1 - | exp 



±V 9 exp 



1 / 2e 3 B M OVT1 



2e 3 B 



V g exp 



-\[^¥ir-R)' 



-^¥(r-R) 



, r<R, 



(20) 



r>R. 



, r<R, 



r>R. 



(21) 



Substituting Eq.(20) into Eq.(ll) in a zero electron mass 
approximation, we get the number density of the electron 



J^Ml-Aexp 



n e (r) = 



2e 3 B 



(r-R) 



,r < R, 



eB 
4 



2e ^(r-R) 



(22) 



r>R. 



Equations (20-22) show that the electrostatic potential, 
electric field, and electron number density are all exponen- 
tial functions of radial distance. They are quite different 
from the classical power-law expressions (jKettner et all 
Il995|) . 

4 Thin crust 

The large outwardly directed electric field is capable of 
supporting some normal material, which then gives birth 
to a thin crust. Two regions should be distinguished just 
outside the quark surface. The first one extends from the 
surface at r = R to the radial distance where the inner nu- 
clear crust (the crust's base) begins, denoted r — R c . The 
associated width, AR = R c — R, is referred to henceforth 
as the gap. The second region begins at R c and extends 
in the radial outward direction toward infinity. Then the 
corresponding Poisson equation reads 



d 2 V 
dr 2 



47re 2 (n e — n q ), r < R, 
= { Aire 2 n e , R < r < R Cl 
4ire 2 (n e - n lon ), R c < r, 



(23) 



where rii on is the positive charge density of ions. For sim- 
plicity and convenience, we assume the electrostatic po- 
tential in the nuc lear crust reg i me (r > R n ) is constant, as 
several authors l|Alcock et all Il986l iKettner et al Il995|) 
have done, i.e., V(r > R c ) = V c (= const). Then the 
boundary conditions become V — > V q as r — > R m and 
V -> V c as r R c . 

Since the mass of ions in the crust is too large for our 
range of magnetic field, the lattice and positive charge 
density of ions should not be influenced by magnetic field. 
For a crust, we define a symbolic electrostatic potential 
V* as 
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Fig. 4. V q , V(R) and V c vs. the magnetic field intensity B. 
A representative value for symbolic electrostatic potential 
V* = lOMeV. 



This new symbol denotes the effective potential at the 
base of the crust in the absence of a magnetic field. We 
use it to designate different crusts. 

The numerical results of V q , V(R), and V c (for V* — 
lOMeV) are shown in Fig. 4. The conventional equation 

V(R) — jV q + jys is also correct for weak magnetic 

fields [B < 10 15 G), while V(R) = \V q + |^ for ultra- 
strong magnetic field approximation. Using this numerical 
result, we show the results of calculations of the electro- 
static potential with respect to different mag netic fields for 
V* = lOMeV in Fig. 5. In previous studies |Alcock et all 
ll986tlKettner et aill995llHuang &Tulll997(l . a minimum 
value of ~ 200fm is established as the lower bound on 
gap width necessary to guarantee the crust against strong 
interaction with the star's strange core. Figure 6 shows 
the curves of gap width versus magnetic field intensity 
for several crusts. If we adopt 200fm as a criterion for 
gap width, then the maximum value of V* may be 9MeV: 
the corresponding crust could be stable for all magnetic 
fields in our range (B < 10 17 G). Our results demon- 
strate that the gap width of various crusts converges to 
a narrow region (~200-400fm) with the increase in B. 
For lighter crusts, the width decreases considerably, i.e. 
from AR = 1556fm at B = 10 15 G to AR = 387fm at 
B = 10 17 G for V* = 3MeV. 

However, a realistic gap width should be influence d 
by the mechanical balance condition l|Huang fc Lulll99^ . 
This condition requires the electron layer above the quark 
surface to partially penetrates the crust. The gap width 
can be over 200fm for only the very lighter crusts when 
the penetrated part is deducted from our present widths. 
The details will be studied in our further work. 



Fig. 5. Same as Fig. 2, but for SSs with a crust with V* = 
lOMeV. 



10 z 




10 15 10 16 10 



B/G 

Fig. 6. Gap width vs. magnetic field intensity. The labels 
refer to symbolic electrostatic potential (in MeV). 

5. Conclusions 

In this paper, the properties of the electric field near 
the surface of strange stars with magnetic field are stud- 
ied thoroughly. We find that a weak (B < 10 15 G) mag- 
netic field has hardly any effect on the surface electric 
field of bare SSs. Secondly, in the middle range (B ~ 
10 15 - 2.5 x 10 16 G) of magnetic field intensity, both the 
electrostatic potential and electric field have an oscillation 
as B changes. Thirdly, for ultra-strong (B > 2.5 x 10 16 G) 
magnetic fields, the electrostatic potential at quark sur- 
face reduces considerably, accompanied by an increase in 
the electric field. The electrostatic potential, the electric 
field, and the electron number density are all exponential 
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functions of radial distance, which is quite different from 
the field-free case. 

We took the effect of magnetic field on the nuclear 
crust into account. With the increase in magnetic field 
intensity, the gap width of various crusts converges to a 
narrow region (~ 200 — 400fm) . The lighter crusts have a 
quick decrease in gap width, while the heavier ones suffer 
a slight increase in gap as B is strengthened. We infer that 
those may influence the mass of the crust. 
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